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APPLICATION  OF  THE  HIGH  TEMPERATURE  CALORIMETER 
TO  THE  DETERMINATION  OF  THE  HEATS  OF  FORMATION  OF  Na-Sn  AND  Li-Sn  ALLOYS 


By  George  R.  Barber,  Leo  Brewer,  LeRoy  A.  Bromley,  and  Raleigh  L.  McKisson 


ABSTRACT 

A  method  has  been  devised  for  using  the  high  frequency  induction  coil  as  a  high  temperature  cal¬ 
orimeter.  The  method  enables  one  to  determine  the  heat  effects  produced  at  high  temperatures.  The 
method  has  been  applied  to  the  determination  of  the  heats  of  formation  of  the  alloys  Na-Sn  and  Li-Sn. 
The  average  heat  of  formation  (AH)  determined  for  Na-Sn  is  -“9.63  ±  0.40  kilocalories  per  mole  at  STS'" 
K.  This  value  becomes  -10.6  t  0.6  kilocalories  per  mole  when  extrapolated  to  298®  K.  The  average 
heat  of  formation  (AH)  determined  for  Li-Sn  is  ~17.57  t  0.83  kilocalories  per  mole  at  850®  K.  This  val¬ 
ue  becomes  -18.62  t  1.50  kilocalories  per  mole  when  extrapolated  to  298®  K.  These  results  are  com¬ 
pared  with  the  results  of  other  investigators. 


PART  I 


INTRODUCTION  TO  THE  WORK  ON  Na-Sn 

The  high  frequency  induction  coil  was  first  used  as  a  calorimeter  by  Chipman  and  Grant.^  Their 
equipment  consisted  of  a  large  silica  crucible  containing  the  charge,  which  was  heated  in  an  induction 
coil.  The  cover  of  the  crucible  was  provided  with  a  small  hole  so  that  additions  and  sampling  could 
be  readily  performed.  An  inert  atmosphere  was  maintained  by  bleeding  argon  into  the  cavity  and  letting 
it  sweep  out  the  hole.  Iron  was  placed  in  the  crucible  and  heated  to  the  desired  temperature,  then 
weighed  samples  of  silicon  and  iron  were  dropped  into  the  melt.  The  iron  was  used  to  estimate  the  heat 
capacity.  The  heat  effect  was  calculated  from  the  heat  capacity  and  the  temperature  difference. 

The  system  that  Chipman  and  Grant  selected  was  an  excellent  choice,  as  the  temperatures  involved 
were  quite  high,  and  considerable  electromagnetic  stirring  was  present.  This  stirring  made  for  rapid 
approach  to  thermal  equilibrium.  The  large  melt  with  its  heat  capacity  prevented  a  great  temperature 
rise,  since  a  relatively  small  sample  was  added.*  Thus  the  temperature  remained  in  a  small  range, 
and  the  heat  losses  could  be  considered  quite  constant  during  a  run.  The  corrections  for  heat  transfer 
were  quite  small,  and  even  though  they  could  not  be  calculated  with  great  accuracy,  the  effect  on  the 
result  was  small. 


♦The  melts  consisted  of  from  65  to  90  poimds  of  molten  iron,  into  which  one  pound  or  less  of 
silicon  was  dropped  per  addition. 
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The  apparatus  used  in  the  present  research  is  similar  to  that  of  Chipman  and  Grant  in  that  high 
frequency  induction  heating  is  used.  However  the  melts  are  much  smaller,*  consequentijrthe  temper¬ 
ature  rises  are  greater,  and  heat  losses  cannot  be  considered  constant. 

The  primary  purpose  of  this  research  was  to  devise  a  method  of  analysis  of  the  temperature -time 
curves  obtained  from  the  apparatus.  The  secondary  purpose  of  this  research  was  to  measure  the  heat 
of  formation  of  Na-Sn.  The  measurement  was  to  be  performed  by  adding  weighed  sodium  pellets  to  a 
molten  tin  melt,  obtaining  a  temperature -time  curve,  and  analyzing  this  curve  by  the  method  devised. 

This  heat  of  formation  has  been  measured  by  several  investigators.  Biltz  and  Holverscheit?  in 
1924  reported  the  heat  of  formation  of  Na-Sn  to  be  eleven  kilocalories  per  mole.  This  result  was  ob¬ 
tained  by  a  differential  solution  method,  using  ferric  chloride  in  about  6N  hydrochloric  acid  as  the  sol¬ 
vent.  In  1928,  Biltz  and  Meyer’*  offered  recalculated  values  of  the  heat  of  formation,  using  what  they 
considered  to  be  a  better  heat  value  for  one  of  the  intermediate  reactions,  and  quoted  a  value  of  six¬ 
teen  kilocalories  per  mole,  Biltz^  again  gives  this  value  in  1937. 

Kubaschewski,®  and  Kubaschewski  and  Seith,®  report  twelve  kilocalories  per  mole  for  this  reaction. 
Their  measurement  was  direct,  the  alkali  metal  being  melted  and  poured  into  the  molten  tin  in  an  argon 
atmosphere,  and  the  mixture  being  quickly  placed  in  a  large  copper-block  calorimeter.  They  calculated 
their  result  from  the  temperature  rise  of  the  calorimeter  and  its  heat  capacity. 

Many  investigators  have  measured  heats  of  formation  of  various  intermetallic  compounds,  and  oth¬ 
ers  have  derived  useful  facts  from  theoretical  considerations.  A  few  references  to  work  of  this  type 
are  listed  in  the  bibliography,  items  10  to  24. 

DESCRIPTION  OF  APPARATUS  AND  OPERATING  PROCEDURE 

The  apparatus  consists  of  a  20  kilowatt  spark-gap  convertert  and  a  calorimeter  which  is  surrounded 
by  one  of  four  inductance  coils  (the  heating  coU).  The  other  three  coils  are  used  to  adjust  the  power  dis¬ 
sipated  in  the  heating  coil. 

The  calorimeter  (see  Figure  1)  consists  of  a  cylindrical  pyrex  tube  which  is  surrounded  by  a  water 
jacket;  a  water  cooled  copper  base  which  supports  the  zircon  stand  upon  which  the  crucible  rests;  and 
a  ‘ ' head which  supports  the  various  devices  which  pass  through  packing  glands  into  the  interior.  These 
devices  are  the  thermocouples,  the  dispenser,  the  dispenser -cooler,  the  operating  rod,  and  the  flexible 
drive  for  the  stirrer.  Two  thermocouples  are  used;  one,  sheathed  in  a  molybdenum  tube,  dips  into  the 
melt;  the  other  passes  through  a  copper  well  to  a  position  near  the  bottom  of  the  dispenser. 

The  dispenser  (see  Figure  2)  was  designed  to  be  an  air-tight  compartmented  container  which  could 
be  loaded  in  a  dry -box  and  transferred  to  the  calorimeter,  thus  preventing,  to  a  great  extent,  any  con¬ 
tact  of  the  sodium  with  oxygen  or  water  vapor.  The  division  of  the  dispenser  into  four  compartments  is 
effected  by  three  pivoted  plates  which  extend  across  the  opening  and  are  held  in  position  by  cams  on  the 
thermocouple  well.  These  cams  are  arranged  so  that  rotating  the  thermocouple  well  in  increments  of 
90*^  causes  a  section  to  dump  its  contents  onto  the  trap  door  forming  the  bottom  of  the  dispenser.  This 
trap  door,  which  is  shielded  from  the  radiation  of  the  hot  melt  by  three  molybdenum  plates  attached  to 
the  bottom,  is  actuated  by  a  rod  passing  through  the  dispenser  cavity.  The  rod  then  passes  through  a 
packing  gland  in  the  top,  and  ends  in  an  eye.  Between  the  top  of  the  dispenser  and  the  rod-eye,  there  is 
a  spring  imder  tension  which  keeps  the  trap  door  closed  tightly.  When  it  is  desired  to  open  the  door,  the 
operating  rod,  whose  end  is  hooked,  can  be  inserted  in  the  eye  and  a  downward  pressure  applied. 


*The  melts  in  this  research  were  about  100  grams  of  Sn,  to  which  from  4  to  10  grams  of  Na  were 
added  per  addition. 

tAjax  20  kw.  Converter,  Ajax  Electrothermic  Corporation,  Trenton,  New  Jersey. 
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Figure  1.  Calorimeter 
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Figure  2.  Air-tight  dispenser  Figure  3.  Dispenser  Cooler 


The  operating  rod  is  also  used  to  rotate  the  thermocouple  well  by  means  of  an  eye  in  the  top  of  the  well. 

The  dispenser  is  enclosed  and  supported  by  the  dispenser -cooler  (see  Figure  3).  This  consists 
of  two  concentric  copper  tubes  with  the  annular  space  closed  at  each  end  and  supported  by  a  copper 
tube  leading  through  the  head.  Concentric  to  and  inside  of  this  tube  is  a  second  tube  which  passes  into 
the  annular  space  to  the  bottom  and  serves  as  the  water  inlet.  Thus  cold  water  is  introduced  at  the  bot¬ 
tom,  swirls  up  to  the  exit  tube  at  the  top,  and  thence  up  inside  the  supporting  tube  to  a  drain.  The  dis¬ 
penser-cooler  can  be  raised  or  lowered  by  sliding  through  the  packing  gland  so  that  when  an  actual 
addition  is  to  be  made  the  dispenser  will  be  near  the  crucible,  and  after  the  addition  is  completed  it 
can  be  raised  away  from  the  hot  crucible  and  thereby  kept  cooler. 

The  stirring  device  (see  Figure  1,  items  3,  12,  20,  24) Consists  of  a  28-volt,  DC,  250-rpm  motor 
coupled  by  means  of  a  sleeved  flexible  shaft  to  the  stirrer .  The  lower  end  of  the  sleeve  is  clamped 
securely  to  the  thermocouple  support  and  provides  a  sort  of  bearing  for  the  shaft.  The  flexible  shaft 
is  welded  to  a  molybdenum  rod  which  is  guided  by  a  molybdenum  sleeve  beating  supported  by  the  ther¬ 
mocouple  well.  The  end  of  the  rod  is  slotted  and  fitted  with  the  stirrer  blades. 

The  crucible  is  surrounded  with  two  0.001 -inch  thick  molybdenum  shields  to  minimize  radiation 
losses.  This  thickness  of  molybdenum  heats  but  very  little  in  the  high  frequency  field*  and  the  outside 
shield  never  reaches  dull  red  heat  even  though  the  inside  of  the  crucible  be  at  bright  red  heat. 


♦The  frequency  is  aboiit  20,000  cycles  per  second. 
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Any  desired  atmosphere  may  be  used  in  the  calorimeter  by  using  tanked  gas.  In  addition,  the  pip¬ 
ing  system  permits  flow  of  gas  into  either  the  top  or  bottom  of  the  calorimeter.  Any  desired  pressure 
may  be  used,  as  the  gas  exit  line  has  a  solenoid  valve  which  is  operated  by  a  mercury -level  contact 
device.  By  adjusting  this  contact  device,  the  pressure  at  which  the  solenoid  valve  releases  is  selected. 
This  system  cannot  be  used  at  present  when  the  stirrer  is  used,  as  there  is  too  great  a  gas  leak  be¬ 
tween  the  flexible  shaft  and  its  sleeve.  In  this  case  the  valve  at  the  tank  is  set  so  that  the  gas  flow  out 
of  the  shaft  is  sufficient  to  prevent  the  diffusion  of  any  air  into  the  apparatus. 

The  electrical  system  is  indicated  in  Figure  4.  The  large  variable  inductance  consists  of  a  sole¬ 
noid  2  feet  in  diameter,  made  of  copper  tubing,  with  shorting  switches  arranged  as  indicated.  Power 


20  KW  SPARK-GAP 
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o- 

o- 
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CONTROLLED 
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MERCURY- 

LEVEL 
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ADJUSTER 


Figure  4,  Schematic  diagram  of  high-frequency  circuit,  and  constant  temperature  mechanism. 
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dissipation  in  the  heating  coil  can  be  increased  by  shorting  out  turns  of  this  coil,  thus  decreasing  the 
total  inductance  of  the  circuit. 

The  heating  coil  consists  of  a  2 8 -turn  solenoid,  6  1/2  inches  in  diameter.  It  is  arranged  on  a  sup¬ 
port  so  that  the  coil  may  be  raised  or  lowered  to  center  the  crucible  in  the  coil. 

The  other  two  coils  make  up  a  temperature  control  device.  The  smaller  coil,  20  turns  on  a  3  1/2 
inch  diameter,  surrounds  a  water  jacketed  iron  vessel  in  which  there  are  a  few  milliliters  of  water. 

The  only  opening  in  the  vessel  is  in  the  bottom,  at  which  point  a  tube  filled  with  mercury  leads  to  a 
pyrex  tube  fitted  with  two  electrodes.  The  pressure  above  the  mercury  in  this  tube  can  be  controlled 
by  a  mercury  bulb  and  reservoir  as  shown  on  the  diagram.  The  two  electrodes  are  part  of  an  electrical 
control  system  which  actuates  a  reversible  motor.  The  motor,  controlled  in  part  by  a  time-delay  re¬ 
lay,  is  coupled  to  the  controlled  variable  inductance  coil  (20  turns  on  an  8 -inch  diameter)  such  that  the 
length  of  the  coil  is  changed  to  respond  to  a  change  in  pressure  in  the  iron  vessel. 

With  this  arrangement  of  coils  and  contacts,  a  small  increase  of  power  increases  the  heat  gener¬ 
ated  in  the  iron  vessel  and  consequently  increases  its  temperature.  This  increase  in  temperature 
causes  an  increase  in  the  vapor  pressure  of  the  water  in  the  iron  vessel  and  raises  the  mercury  in 
the  contact  chamber.  When  the  contact  circuit  is  completed,  the  motor,  which  normally  decreases  the 
inductance  of  the  controlled  variable  inductance  coil,  reverses  and  increases  its  inductance.  This  in¬ 
creases  slightly  the  power  dissipated  in  this  inductance  and  consequently  reduces  the  power  dissipated 
in  the  other  coils  in  the  system.  This  decrease  in  power  lowers  the  temperature  and  vapor  pressure 
of  the  water  and  ultimately  causes  the  mercury  circuit  to  break.  At  this  time  the  motor  reverses,  and 
the  controlled  variable  inductance  decreases.  This  decrease  results  in  an  increase  in  power,  increases 
the  vapor  pressure  of  the  water  in  the  iron  vessel,  and  closes  the  mercury  contact.  Then  the  cycle  re¬ 
peats.  The  entire  cycle  causes  a  temperature  variation  in  the  melt  of  about  1°C. 

The  time -delay  relay  on  the  motor  was  necessary  because  the  system's  reaction  time  was  too  slow. 
If  the  motor  ran  continuously  it  would  cause  the  coil  to  over -correct  and  no  control  would  be  effected. 
Satisfactory  operation  is  obtained  with  a  relay  setting  to  give  a  15 -second  '^on"  period  and  a  5 -second 
"'off"  period. 

Thus,  this  system  controls  the  power  dissipation  in  the  entire  high  frequency  circuit,  and  various 
levels  of  power  can  be  maintained  merely  by  adjusting  the  height  of  the  mercury  bulb. 

To  operate  the  apparatus,  one  loads  the  dispenser  with  weighed  quantities  of  sodium,  assembles 
the  calorimeter,  and  starts  the  converter.  When  the  tin  has  melted,  the  stirrer  and  the  thermocouple 
are  set  in  place.  Temperatures  are  recorded  on  the  micromax,*  and  when  the  system  is  near  the  de¬ 
sired  temperature,  the  power  setting  and  incremental  inductance  are  adjusted  so  that  the  system  ap¬ 
proaches  equilibrium.  The  control  device  is  then  set  in  operation.  When  the  temperature  of  the  melt 
remains  constant  for  a  few  minutes,  the  temperature  of  the  dispenser  is  taken  and  the  sodium  addition 
made.  Before  and  as  soon  as  possible  after  the  addition,  millivolt  readings  are  taken  on  the  potentiom- 
etert  to  check  the  values  recorded  by  the  micromax. 

When  the  temperature  has  reached  its  new  equilibrium,  the  power  may  be  adjusted  to  yield  the  de¬ 
sired  temperature  for  the  next  addition.  This  procedure  is  repeated  for  all  additions.  When  the  run  is 
completed,  the  calorimeter  is  opened  and  the  dispenser  examined  to  see  if  all  the  pellets  were  added. 
The  melt  is  removed,  ground  in  the  dry  box,  and  anaylzed.  From  the  analysis,  the  amount  of  sodium 
vaporized  is  determined. 


*Micromax,  Leeds  and  Northrup  Company,  Philadelphia. 
tPrecision  potentiometer,  Leeds  and  Northrup  Company,  Philadelphia. 
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DERr/ATION  OF  NEW  METHOD  OF  ANALYSIS  OF  TEMPERATURE -TIME  CURVES 

The  principal  data  received  from  a  run  consist  of  a  temperature -time  curve  traced  by  the  micro¬ 
max,  and  in  addition,  a  few  points  taken  on  the  potentiometer  as  checks.  It  is  therefore  necessary  to 
devise  a  method  of  interpretation  of  these  curves.  Typical  curves  are  shown  in  Figure  5. 

The  rate  of  heat  loss  at  temperature,  t,  can  be  expressed  by  the  following  equation*: 

-qit=aAie(t4-T^) +A2h(t-Ti.),  (1) 

in  which  the  coefficient  for  heat  transfer  by  conduction  is  combined  with  that  for  convection.  When  the 
system  reaches  a  steady-state  temperature,  t,  the  thermal  energy  lost  to  the  surroundings  is  equal  to 
the  thermal  energy  gained  from  the  high  frequency  heating  coil. 

The  rate  of  heat  loss  or  gain  at  any  temperature,  T,  is, 

-qiT=<^Ai€(T4-T4)  +A2h(T-Tr).  (2) 

The  difference  in  rates  of  heat  loss  at  T  and  t  is  the  net  rate  of  heat  transfer  at  T,  and  is  represented 

-  (q^T-  qit)  =  (T^  - 1^)  +  AgCT  - 1).  (3) 

(t4  expanded  in  powers  of  (T  -t),  thus: 

(t4  ^ =  4t3(x  - 1)  +  6t2(T  - 1)^  +  4t(T  ~  t)^  +  (T  - 1)^.  (4) 


Figure  5.  Typical  temperature -time  curves. 


♦Definitions  of  symbols  following  end  of  Part  n. 
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For  the  case  in  which  t  equals  873°K.,  and  (T  -t)  is  125°K.,  the  contributions  to  the  total  of  the  vari¬ 
ous  terms  in  the  expansion  are:  80.9%,  17.4%,  1.66%,  and  0.04%.  For  the  case  in  which  t  equals  873° 
K.,  and  (T  -t)  is  50°K.,  the  contributions  to  the  total  of  the  various  terms  in  the  expansion  are:  91.8%, 
7.9%,  0.3%,  and  less  than  0.1%.  From  the  above,  it  is  evident  that  if  one  limits  (T  -t)  to  less  than 
125° K,  neglecting  the  last  two  terms  in  the  expansion  will  introduce  less  than  2%  error.  Therefore,  an 
allowable  substitution  for  (T^  -t^)  in  equation  (3)  is  4t^(T  —t)  +  6t^(T  —  t)^.  Rewriting  equation  (3), 

-(qiT  -  qit)  [4t3(T  -t)  +  6t2(T  -t)2]  +  AgM?  -t).  (5) 

The  corresponding  equation  written  for  any  temperature,  T,  and  the  final  equilibrium  tempera¬ 
ture,  t',  is: 

-(qiT  f  -t')  +  6t '2(T -t')2]  +A2h(T-t').  (6) 

Letting  aA^e  =  B,  A2h  =  D?  AT '  =  (T  - 1' ),  equation  (6)  becomes,  upon  re¬ 

arranging, 

f  =  (4Bt'^  +  D)AT'  +  6Bt'2(AT')^.  (7) 

Now  let  (4Bt'2  +  D)  =  P,  and  let  6Bt'2  =  R,  then, 

f  =  PAT'  +  R(AT')2.  (8) 

K  we  can  evaluate  P  and  R,  then  f  can  be  calculated  for  a  number  of  (AT')'s  and  an  f  curve  plotted 
from  which  one  can  read  values  for  the  rate  of  heat  transfer  for  the  system  at  any  temperature  differ - 


The  change  in  convection  characteristics  and  the  change  in  radiation  characteristics  are  consider¬ 
ed  small  for  any  run,  therefore  values  of  P,  or  P',  and  R,  or  R',  calculated  for  a  portion  of  the  run,  are 
representative  of  the  entire  run.  Thus  f  values  relative  to  any  steady -state  temperature  may  be  read 
from  the  f  curve  mentioned  above  when  temperature  differences  are  used  as  the  basis  of  comparison.* 

By  plotting  f  (cal/min)  against  9  (min),  one  obtains  areas  directly  representing  quantities  of  heat. 
This  fact  will  be  used  later  in  the  discussion. 

To  evaluate  P  and  R,  one  can  use  the  cooling  portion  of  the  temperature -time  curve  beyond  point 
'‘a”.  Figure  6,  if  one  assumes  this  curve  to  approach  a  theoretical  cooling  curve  type.  Then, 

f  =  PAT'  +  R(AT')2;  (8) 

dQ=-£de;  (9) 

dQ  =  £  wc„  dT  =  C  dT;  (10) 

(11) 

de  de  ’ 


and  ^ 
dO 


[PaT'  +R(aT')2]  . 


Rearranging  for  integration, 


-T  dT'  _  rdS  . 

J  R(AT')2 +PAT'  J  Cp 


♦For  the  case  of  t;|^  =  873°K  and  t2  =  890°K,  a  comparison  of  Fs  for  equal  values  of  B,  D,  and  T', 
is:  for  tj,  214;  for  t2,  220.  The  difference  is  less  than  3%.  B  was  taken  as  2(10"^),  D  as  3.0,  and  the 
temperature  difference  as  25° K. 
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Since  AT'=  (T  —  t'),  d(AT')  =  dT.  Substituting  and  integrating, 


P  +  RAT ' 
AT' 


Then,  solving  for  6,  and  rearranging, 

S,  _ 

0  =  —  p  1  +  5.AT'  P  +  C, 
P 


(14) 


(15) 


or, 


1  +^AT' 


+  c. 


(16) 


If  the  region  of  the  cooling  curve  does  approach  a  theoretical  cooling  curve,  one  obtains  a  straight 
line  whose  slope  is  “P/Cp  by  plotting 


In - - - against  9, 

1 


Since  R  and  P  are  both  unknown,  this  plot  cannot  be  made  directly.  However,  as  a  first  approximation, 
one  can  plot  In  AT'  vs.  9  and  obtain  a  value  from  the  slope,  from  which  P  can  be  calculated  when  Cp 
is  known.  This  value  should  be  quite  close  to  the  correct  value  of  P,  if  R  is  small  relative  to  P*,  as 
can  be  seen  in  equation  (16). 

Using  this  value  of  P,  one  can  solve  for  B  from  the  definition  P  =  4B(t')^  +  D,  knowing  D  (Appendix 
I).  With  B,  R  can  be  calculated  from  the  relation  R  =  6B(t')^.  These  approximate  values  of  P  and  R 
now  permit  calculation  and  plotting  of 

AT ' 

In - - — -vs.  0. 

l+fAT' 

From  the  slope  of  the  smoothed  curve,  new  values  P '  and  R'  can  be  calculated.  This  procedure  may 
be  repeated  until  sufficient  accuracy  is  obtained.  The  values  of  P  and  R,  necessary  for  plotting  the  f 
curve,  are  now  found. 

It  is  desirable  to  be  able  to  extrapolate  the  cooling  curve  obtained  from  the  experimental  data  to 
such  time  that  the  AT  can  be  evaluated  as  if  the  reaction  had  proceeded  instantaneously.  This  is  most 
easily  accomplished  by  making  use  of  the  final 

^T' 

In - - - vs.  0. 

1  +  -AT ' 


curve  obtained  above.  Select  values  of 


In 


AT'' 

1 +I7AT'' 


from  the  straight  line  determined  by  values  of  the  similar  expression  in  AT'.  Solve  for  values  of  aT' 
and  plot  the  corresponding  extrapolated  T’s  on  the  T  vs.  0  graph,  as  indicated  on  Figure  6.  K  one  can 
obtain  the  0j  corresponding  to  an  instantaneous  reaction,  then  one  can  also  obtain  the  accompanying 
aT  from  the  extrapolated  curve  above. 


*The  calculations  show  P  to  be  about  lo,  and  R  to  be  about  10 “3. 


Ill 


Figure  6.  Temperature -time  curve,  Run  1,  Curve  2,  Na 
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Superposing  an  ideal,  instantaneous  T  vs.  0  curve  on  a  typical  T  vs.  0  curve  defines  areas  which 
represent  relative  heat  gains  and  losses  by  the  system.  (Areas  I,  II,  and  in  on  Figure  6.)  Those  areas 
in  which  the  ideal  curve  lies  above  the  actual  are  '‘heat  gain’’  areas,  since  the  system  is  losing  heat 
at  a  lower  rate  than  it  gains  heat  (T<t,  area  I);  or  the  system,  although  losing  heat  at  a  higher  rate  than 
it  gains  heat,  is  losing  heat  at  a  lower  rate  than  the  ideal  system  would,  resulting  in  a  net  equivalent 
gain  of  energy  (area  III).  Those  areas  in  which  the  ideal  curve  lies  below  the  actual  curve  are  heat 
loss”  regions,  because  in  this  case  (area  II)  the  system  loses  heat  at  a  higher  rate  than  it  gains  heat, 
resulting  in  a  net  loss  of  energy.  The  proper  0,  0^,  for  the  instantaneous  reaction  is  the  one  for  which 
the  heat  losses  exactly  equal  the  heat  gains.  This  0  is  obtained  by  trial -and~err or  balancing  of  Area  I 
plus  Area  HI  against  Area  II  on  an  f  plot. 

In  the  construction  of  the  f  plot,  an  f  curve  should  be  used  which  was  calculated  for  the  system  in 
a  completely  nonreactive  state,  for  any  other  f,  such  as  one  determined  from  values  of  P'  and  R'  for 
a  sodium  addition,  is  a  function  of  the  rate  of  reaction.*  The  only  f  that  is  not  a  function  of  the  rate 
of  reaction  is  the  one  determined  from  the  molybdenum  addition. 

The  f  plot  is  constructed  by  reading  T’s  for  values  of  0  on  the  T  vs.  0  curve,  calculating  the  cor¬ 
responding  T  —  t,  reading  the  values  for  f  on  the  f  curve,  and  plotting  the  value  of  f,  at  the  proper  0, 
on  the  f  plot  vs.  0  plot.  Figure  7. 

Since  the  units  of  f  are  calories  per  minute,  an  area  produced  by  plotting  f  vs.  0  in  minutes  re¬ 
presents  a  quantity  of  heat.  In  the  ideal  reaction,  there  would  be  no  loss  or  gain  of  heat  by  heat  trans¬ 
fer  because  of  its  instantaneous  character.  It  is  therefore  necessary  to  choose  the  0j  as  the  0  for 
which  the  net  heat  transferred  in  the  actual  case  is  zero.  The  physical  significance  of  the  ideal  plot 
is,  for  convenience  in  graphical  calculation,  to  combine  an  instantaneous  reaction  in  which  equilibrium 
is  partially  reached  with  a  subsequent  reaction  period  during  which  equilibrium  is  quantitatively  ap¬ 
proached. 


♦The  rate  of  heat  transfer,  f,  at  any  time,  0,  is  the  sum  of  two  rates,  thus:  =  fj-eact. 

^heat  trans  *  ^i^^e  we  wish  to  balance  heat  losses  and  gains  due  to  heat  transfer,  we  must  not  use  an 
f  which  may  have  an  f^eact  included,  therefore,  we  must  use  the  f  determined  from  the  character¬ 
istics  of  a  molybdenum  addition. 


Figure  7.  f  curve,  and  f  j)fot. 


12 


AECD  -  2294 


While  Area  I  is  bounded,  so  that  it  represents  a  definite  amount  of  heat  absorbed,  Area  II  is  not 
bounded  on  the  right.  The  aim  is  to  locate  the  boimd  of  Area  n  such  that  the  effective  net  heat  trans¬ 
ferred  from  the  system  (represented  by  Areas  I,  II,  and  III)  is  zero. 

Now  let  us  consider  the  temperature -time  curve  of  Figure  6.  If  the  curve  beyond  point  a  is  a  typ¬ 
ical  cooling  curve -type,  which  can  be  determined  by  plotting  on  a  In  AT'  vs.  0  plot,  then  it  is  also  a 
part  of  the  ideal  temperature  curve,  and  as  such,  can  be  extrapolated  to  the  left.*  Since  we  are  trans¬ 
ferring  the  actual  T  vs.  0  curve  to  an  ideal  curve,  there  is  a  point  on  the  extrapolated  curve  which  re¬ 
presents  the  highest  temperature  that  the  system  would  attain  in  the  ideal  reaction.  This  temperature 
is  defined  by  the  condition  that  net  heat  transfer  at  the  corresponding  time  is  zero. 

The  extrapolated  T  vs.  0  curve  and  the  actual  temperature -time  curve  partially  define  a  third  area. 
Area  in,  which,  when  plotted!  on  an  f  plot,  represents  a  quantity  of  heat  gained  by  the  system.  Area 
in,  like  Area  II,  is  unbounded.  When  these  two  areas  are  plotted  at  their  respective  times,  as  shown 
on  Figure  7,  a  unique  bounding  vertical  line  can  be  found  such  that  the  sum  of  Area  I  plus  Area  in 
equals  Area  II.  This  is  exactly  the  condition  placed  on  the  0^,  and  since  the  line  is  unique,  it  deter¬ 
mines  0j.  Then  the  point  of  intersection  of  the  vertical  0^  line  with  the  extrapolated  temperature  curve 
defines  the  temperature  which  an  ideal  reaction  would  have  attained.  The  temperature  difference  be¬ 
tween  this  derived  temperature  and  the  initial  steady-state  temperature  defines  a  quantity  of  heat,  Qj, 
from  the  relation, 

Ql  =  CpAT.  (17) 

The  Cp  that  should  be  used  in  equation  (17)  is  the  Cp  of  the  system  after  the  sodium  has  been 
added.  For  runs  in  which  multiple  additions  are  made,  values  of  Cp  must  be  calculated  for  all  sodium 
additions  but  the  last.  This  is  done  by  assuming  that  the  Cp  for  the  system  at  any  time  is  equal  to  L 
Cp  less  a  factor,  F,  times  the  Cp  of  all  the  sodium  added  after  the  addition  in  question.  The  evaluation 
of  F  is  presented  in  Appendix  II. 

In  the  event  that  the  reaction  is  not  complete  at  this  point,  and  heat  is  being  liberated  throughout 
the  ^ ‘cooling^’  part  of  the  T  vs.  0  curve,  an  estimation  of  this  added  heat  must  be  made. 

The  criterion  for  determining  whether  the  reaction  was  complete  at  '^a'"  is  whether  the  final  slope 
of  the 


In  — - vs  0 

l+f  AT' 

plot  for  a  sodium  addition  is  less  than  the  corresponding  slope  for  the  molybdenum  addition  at  the  end 
of  the  run.  If  the  slope  is  materially  less,  the  reaction  was  not  complete  at  a.  Any  unreacted  sodium 
present  at  point  a  will  liberate  its  heat  of  reaction  before  the  system  reaches  its  final  equilibrium 
temperature. 

At  this  point  in  the  development,  the  T  vs.  0  curve  obtained  experimentally  has  been  transformed 
to  an  equivalent  theoretical  instantaneous  reaction  T  vs.  0  curve,  and  we  have  accoimted  for  some  frac¬ 
tion  of  the  added  sodium  which  can  be  assumed  to  have  reacted  instantaneously  to  arrive  at  point  C, 
Figure  6.  The  characteristics  of  a  cooling  curve  of  a  completely  reacted  system  can  be  approximated 
from  the  molybdenum  addition,  in  which  addition  no  reactive  material  is  introduced.  By  passing  such 
a  curve  through  point  C  so  that  its  equilibrium  temperature,  at  0  equals  infinity,  is  t'  an  area  is  de¬ 
fined  (Area  IV)  which  represents  the  quantity  of  heat  liberated  by  the  remaining  fraction  of  the  sodium 


♦The  method  can  be  modified  to  apply  where  the  ideal  curve  is  not  approximated. 

tThe  ordinates  of  this  area  are  plotted  as  the  difference  between  the  f  for  the  temperature  on  the 
extrapolated  curve  and  the  f  for  the  temperature  on  the  actual  curve. 
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before  the  system  reaches  t' .  This  quantity  of  heat  may  be  determined  by  plotting  the  area  on  an  f ' 
plot,  on  which  areas  are  proportional  to  heat.  Any  ordinate  of  the  corresponding  area  on  the  f '  curve 
is  found  by  subtracting  the  f  determined  for  temperatures  on  the  lower  (molybdenum -type)  curve  from 
the  f  determined  for  temperatures  on  the  upper  (actual)  curve,  both  temperatures  being  taken  at  the 
same  0. 

Since  it  is  impossible  to  measure  an  area  out  to  infinite  time,  the  area  can  be  broken  up  into  two 
parts  for  actual  measurement.  The  direct  measurement  of  area  on  the  f '  plot  should  be  extended  to 
some  practical  limit  at  which  time  it  may  be  safely  assumed  that  the  reaction  is  complete,  and  that 
the  T  vs.  0  curve  is  an  exact  cooling  curve  beyond  this  limit.  The  heat  represented  by  this  measur¬ 
able  area  is  designated  by  the  symbol  Qjj.  The  heat  represented  by  the  area  out  to  infinite  time  can 
be  calculated  analytically,  thus: 


dQ  =  “f  d0, 

(9) 

and, 

f  =  PAT'  +  R(AT')^. 

(8) 

For  small  (AT'), 

f  =  PAT', 

(18) 

and, 

dQ  =  ■-PAT'd0 

(19) 

From  equation  (12),  when  AT'  is  small, 

dT  ^  ^  n., 

d9  ’ 

(20) 

and. 

P 

d0  "  Cp 

(21) 

Transforming,  and  integrating, 

r  r  p  le 

[inAT'l  =  -— 0 

(^T')o  LCpJ^^ 

(22) 

When  0Q  is  taken  to  be  zero,  and  AT  '  at  0q  equals  (AT')o> 

in  _  _  P  , 

(at')o  Cp® 

(2?) 

and 

(AT')g  =  (AT')(,e-(P/Cp)e- 

(24) 

Let  Qjjj  be  the  heat  in  question,  then. 

dQjii  =  (dQ^ctual  curve  ”  ^^theo.  curve) » 

(25) 

then,  using  equation  (19), 

actual  Qurve  ^’^^Hheo.  curve) 

(26) 

Substituting  equation  (24)  in  equation  (26)  and  integrating, 

I  <lQin  =  -  P 


14 


AECD  -  2294 


+  /P{^T')otheo.  de,  (27) 

0 

or, 

^ni  “  ^)o  actual  "  )q  theoj  ’  (28) 

which,  by  definition,  is  set  equal  to, 

Qm  =  Cp  AT  ^  (29) 

This  simple  expression,  then,  gives  a  means  of  measurement  of  the  area  from  some  point  of  small 
(AT')^s  to  infinite  time. 

To  account  for  the  heat  lost  from  the  system  due  to  vaporization  of  sodium,  one  must  first  deter¬ 
mine  the  amount  of  vaporization  occurring  during  the  particular  addition.  The  total  vaporization  can  be 
calculated  from  an  analysis  of  the  final  melt  since  the  extremely  low  vapor  pressure  of  tin  prevents  its 
vaporization.  The  problem  then  becomes  one  of  determining  the  fraction  of  the  total  sodium  vaporized 
that  was  vaporized  in  the  particular  addition. 

It  is  assumed  that  since  sodium  has  a  much  lower  density  than  tin,  a  phase  of  molten  sodium  tends 
to  cover  the  tin  phase  until  the  sodium  reacts.  Therefore  the  ratio  of  duration  of  a  particular  addition 
to  the  total  for  all  additions  of  a  run  is  used  to  proportion  the  amount  of  sodium  vaporized  during  the 
addition.  The  quantity  of  heat  lost  by  vaporization  is  designated  by  the  symbol  Qq,  and  can  be  calculated 
from  the  known  heat  required  to  heat  solid  sodium  at  the  temperature  of  the  dispenser  to  sodium  vapor 
at  the  temperature  of  the  melt. 

The  total  amount  of  heat  liberated  in  any  addition  is  the  sum  of  the  various  Q’s  as  calculated  above. 
The  AH  for  the  reaction,  Na^^^i^Cs)  +  Snj^Q^(l)  NaSnj^Q|.(l),  is  this  summation  of  Q’s  (XQ)  divided  by 
the  moles  of  sodium  added.  Then  the  AH  of  the  reaction  Na^ot^^)  NaSn^oti^)  can  be  cal¬ 

culated  by  adding  to  the  AH  above  the  heat  required  to  raise  the  sodium  from  the  dispenser  temperature 
to  the  temperature  of  the  melt,  t. 


SAMPLE  CALCULATION 

The  sample  calculation  presented  is  of  the  third  sodium  addition  of  the  first  run.  The  molybdenum 
addition  is  calculated  first,  to  find  the  Cp  of  the  system. 


Table  1.  Data  for  In  AT^  vs.  0  plot,  run  L  curve  4,  Mo. 


j.  clulk:  j.  .  J. 

e 

AWA  All  t 

AT' 

AT'  * 

AT' 

1  +  (R/P)AT' 

9 

2.3 

- 

2.3 

8 

3.2 

- 

3.2 

6 

6.1 

6.1 

6.1 

4 

11.0 

11.0 

10.9 

2 

19.6 

19.9 

19.2 

1.5 

23.3 

- 

22.8 

1 

27.8 

26.8 

27.1 

0.5 

- 

31.0 

- 

0 

- 

36.0 

- 

♦These  values  are  read  from  the  line  drawn  through  the  (AT')’s  to  extrapolate  T’s  to  the  0’s 
indicated. 
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In  the  calculation  of  X  Cp,  the  interpolated  and  extrapolated  values  of  T'  *  above,  which  were  read 
from  the  line  on  the  In  AT'  vs.  6  plot,  Figure  9,  are  plotted  on  the  T  vs.  6  plot,  Figure  8,  and  a  smooth 
curve  is  drawn  through  these  points.  By  balancing  the  areas  on  the  T  vs.  0  plot  as  if  they  were  on  an 
f  plot,  the  approximate  time  0^  and  the  corresponding  AT  can  be  found.  This  value  of  AT  permits  cal¬ 
culation  of  a  X  Cp.  It  was  found  that  6^  was  0.27  minutes  and  AT  was  (614  -  575)  or  39.0°C.  From  heat 
capacity  data  given  by  K.  K.  Kelley,’’^  the  aH  for  molybdenum  is  calculated: 

AH  =  5.69(848  -  318)  +  0.94(10-3)(8482  -  gisS)  +  (5.03)(10‘^)/848 
-  {5.03)(1o4)/318; 

=  3020  +  581  +  59.4  -  158  =  3502.4  or  3.50  kcal/mole 

The  weight  of  molybdenum  was  25.85  grams  or  0.269  moles.  Then  £  Gp  =  (0.269)(3.50)/(39;())  =  0.02414 
or  24.14  cal/  “C.  Using  this  value  for  £  Cp,  P  can  be  calculated;  P  =  (0.291)(24.14)  or  7.03.  Using  the 
value  for  D,  2.62,  as  calculated  in  Appendix  I,  and  the  value  for  P,  B  is  found  from  the  equation,  B  = 
(P-D)/4(t')3  =  (7.03  -2.62)/4(881)3,  or  B  =  1.61(10‘9).  Then,  since  R  =  6B(t')2,  R  =  6(1.61)(l0'9)(88l)2 
or  R  =  0.0075.. 


These  values  of  P  and  R  allow  the  calculation  of  the 

AT' 

1+|at' 

column  of  Table  1.  These  points  are  plotted  as  shown  bn  Figure  9.  The  Slope  of  the  line  does  not 
change j  and  the  values  of  P  and  R  were  used  to  calculate  an  f  curve  and  make  an  f  plot,  but  since  the 
was  not  changed  from  its  original  value,  these  graphs  are  not  shown.  This  effect  is  due  to  the  fact  that 
the  areas  are  small,  and  the  (AT'  )’s  are  small. 

Calculation  of  Sodium  Vaporized  in  Run  1 

The  molal  ratio  of  tin  to  sodium  from  amounts  added  for  the  run  is  equal  to  (22.997)(95.8)/(118.7) 
(17.97)  or  1.033.  The  molal  ratio  of  tin  to  sodium  from  analysis  is  1.060.  The  amount  of  sodium  in  the 
final  melt,  assuming  that  no  tin  is  vaporized  during  the  run,  is  (1.033)(17.97)/(1.060)  or  17.51  grams. 
Therefore  the  sodium  vaporized  is  (17.97  -  17.51)  or  0.46  g.  This  total  amount  of  vaporization  is  pro¬ 
rated  in  proportion  to  the  time  of  an  addition,  thus: 

Time,  Curve  1  (20.0  -0.8)  =  19.2  min 
Curve  2  (23.1  - 0.6)  =  22.5  min 
Curve  3.  (25.0-0.1)  =  24.9  min 
Total  66.6  min 

Amounts  vaporized  in  each  addition  are: 

Curve  1,  (19.2)(0.46)/(66.6)  =  0.133  g 
Curve  2,  (22.5)(0.46)/(66.6)  =  0.155  g 
Curve  3,  (24.9)(0.46)/(66.6)  =  0.172  g 

Actual  amOimts  of  sodium  added  per  addition: 

Curve  1,  (8.27  -  0.133)  =  8.137  g 
Curve  2,  (5.20  -0.155)  =  5.045  g 
Curve  3,  (4.50  -  0.172)  =  4.328  g 


♦These  values  are  read  from  the  line  drawn  through  the  (AT')'s  to  extrapolate  T's  to  the  0's 
indicated. 
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Figure  8.  Temperature -time  curve,  Run  1,  Curve  4,  Mo. 
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Calculation  of  the  Heat  of  Reaction  from  Curve  3,  Run  1 

The  weight  of  sodium  added  in  this  addition  is  4.328  grams.  The  Cp  for  the  system  for  this  addition 
is  24.14  cal/  **C.  The  equilibrium  temperature  for  this  run  was  found  to  be  613‘'C,  t' ,  and  D  was  2.62. 

The  values  of  data  given  in  Table  2  for  this  addition  give  a  slope  on  the  In  AT '  vs.  0  plot  of  “0.108. 
Using  this  slope,  and  the  value  for  the  Cp  of  the  system,  P  is  (0.108) (2 4. 14)  or  2.61.  Then  B  is  (P  ~  D)/ 
4(t'  r,  or  when  the  values  are  substituted,  (2.61  “  2.62)/4(886)3  or  “3.60(10“^^).  R  is  6B(t')^  =  6(~3.60) 
(10'^2)(886)2  or  -1.10(10“^)  * 


Table  2.  Data  for  In  AT^  vs.  $  plot,  run  1,  curve  3,  Na. 

AT^  AT^  ^ _ 

AT' 


e 


1  +|r^" 


AT' 


27 

0.9 

0.9 

- 

- 

23 

2.1 

2.1 

- 

- 

20 

3.4 

3.4 

- 

- 

17 

5.3 

5.3 

- 

- 

15.5 

6.4 

6.4 

- 

- 

14 

7.6 

7.6 

7.5 

7.5 

12.5 

8.8 

8.8 

8.8 

8.8 

11 

10.1 

10.0 

10.3 

10.3 

10 

- 

- 

11.6 

11.6 

9 

- 

12.8 

12.8 

Since  R  is  very  small,  the  values  of  , 

1+|aT' 

are  practically  the  same  as  the  corresponding  values  of  AT' .  The  slope  remains  the  same,  and  P  =  P' , 
B  =  B',  and  R  =  R'.  Values  of 

1 +I7AT" 

read  from  the  plot  and  the  corresponding  (AT")'s,  are  tabulated  in  Table  2  above.  The  T's  calculated 
from  the  (AT '  ')'s  are  plotted  on  the  T  vs.  6  plot,  Figure  10,  An  f  curve  must  now  be  plotted  from  the 
relation,  f  =  P'A  T'  +  R'  (AT'  )^,  where  P'  and  R'  are  the  values  from  the  molybdenum  addition.  Data 
appear  in  Table  3.  The  f  plot  is  constructed  using  values  of  (T  “t)  read  from  the  T  vs.  0  curve  to  find 
the  correct  value  of  f  from  the  f  curve,  and  then  plotting  the  values  of  f  against  the  proper  0.  This  plot 
appears  in  Figure  13.  Areas  are  balanced  on  this  f  plot  as  indicated,  and  the  0j^  is  found.  At  this  0^, 
the  AT  is  read  from  the  T  vs.  0  curve,  for  the  ideal  reaction,  and  is  11.6°C.  at  0|  =  9.85  minutes 
(Figure  10). 


*It  is  possible  for  R  and  B  to  be  negative  for  a  sodium  addition  when  there  is  a  great  deal  of  so¬ 
dium  reacting  in  the  latter  part  of  the  addition.  The  negative  value  occurs  because  the  slope  of  the  T  vs. 
0  curve  is  quite  flat  in  cases  of  this  type. 
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Figure  10,  Temperature -time  curve.  Run  1,  Curve  3,  Na. 


Figure  11.  Ln  AT'  vs.  6  plot,  Run  1,  Curve  3,  Na. 
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X  u. 

AT' 

f 

AT' 

f 

50 

370.0 

0 

0 

30 

217.5 

-15 

-104.1 

15 

107.1 

-30 

-204.5 

The  value  of  Qj  can  now  be  calculated.  Qj  is  equal  to  Cp(AT)  =  (24.14) (11.6)  or  280  calories. 

Using  data  from  K.  K.  Kelley,’  Qq  can  be  calculated. 

Qq  =  (0.172/22.997)  [5.01(371  -  318)  +  (5.36/2)(10’3)(37l2  -  318^)  + 

630  +  7.50(886  -  371)  +  23890®] 

=  (0.0075)(266  +  98  +  630  +  3863  +  23890)  or  0.216  kcal. 

The  calculation  of  Qjj  involves  plotting  an  f '  plot.  This  involves  plotting  f's  from  the  table  above, 
as  determined  from  the  (at')'s  of  area  IV,  Figure  10,  against  the  corresponding  0’s,  The  resulting 
f '  plot  is  shown  in  Figure  14.  The  area  was  measured  and  found  to  be  1004  units,  at  5/16  cal/unit, 
which  equals  314  calories. 

The  calculation  of  Qjjj  consists  of  calculating  the  quantity  (AT)  from  the  T  vs.  0  plot,  and  sub¬ 
stituting  in  the  equation  Qjjj  =  ^p(AT)  ;  therefore  equal  to  (24.14) (1.2)  or  29  calories. 

Qo  +  Qi  +  Qii  +  Qin  “  (0*216  +  0.280  +  0.314  +  0.029)  or  0.839  kilocalories  for  4,367  grams  of 
sodium. 

(Qq  +  Qi  +  Qn  +  Qni)  “  (0.839) (22. 997)/(4.328)  or  4.458  kilocalories  per  mole. 

Q  (266  +  98  +  3863  +  630)  or  4.857  kilocalories  per  mole. 

-  XQ  is  equal  to  the  A  H  for  the  reaction, 

Nagge  +  Sngse  NaSn886» 

and  this  AH  is  -<4.458  +  4.857)  or  -9.315  kilocalories  per  mole. 

For  the  reaction, 

Na873  +  Sn873  —  NaSn873, 
the  A  H  from  the  third  sodium  addition  is: 

-9.315  -  (1  -0.8)(7.5)(886  -  873)(10“^),  or  -9.33  kilocalories  per  mole. 

TABULATION  OF  RESULTS  ON  Na-Sn 

In  calculating  the  final  result,  the  three  additions  of  run  1  were  averaged,  the  three  additions  of 
run  10  and  run  11  combined  were  averaged,  and  the  average  of  these  two  averages  taken  as  the  result. 
This  is  indicated  on  Table  5. 

Errors  were  calculated*  for  each  type  of  run,  and  it  was  assumed  that  the  same  error  was  pre¬ 
sent  in  all  runs  of  the  same  type. 


♦Errors  were  calculated  on  run  1,  addition  2;  and  on  run  11,  addition  1. 


f  IN  CAL/MINUTE  CAL/MINUTE 
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Figure  13.  f  plot,  Run  1,  Curve  3,  Na. 


9  12  15  18  21  24  27 

6  IN  MINUTES 

Figure  14.  f'  plot,  Run  1,  Curve  3,  Na. 
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Table  4,  Summary  of  runs  and  results. 


Run 

No.  Na 
additions 

Heat  of 
reaction 
kcal/ mole 

Comments 

1 

3 

1st  -9.77 

No  stirrer. 

2nd  -  10.76 

No  stirrer. 

3rd -9.33 

No  stirrer. 

2 

3 

Dispenser  failed. 

3 

3 

- 

Vaporization  losses  excessively  high. 

4 

- 

Dispenser  failed. 

5 

3 

- 

Third  addition  incomplete. 

6 

- 

Deflecting  shield  failed. 

7 

- 

- 

Crucible  leaked. 

8 

2 

~ 

Vaporization  losses  excessively  high. 

9 

2 

- 

Deflecting  shield  failed. 

10 

1 

-8.87 

With  stirrer. 

11 

2 

1st -9.73 

With  stirrer. 

2nd  —  9,30 

With  stirrer. 

Table  5.  Results  of  determinations,  Na-Sn. 

Rim 

Addition 

Result 

Kcal/mole  at  873°  K 

1 

1 

-9.77  i  0.77 

2 

-10.76  t  0.77 

3 

-9.38  t  0.77 

Average 

-9.95  +0.45 

10 

1 

-e.87  ±  1.14 

11 

1 

-9.73  +  1.14 

2  -«.30  +  1.14 

Average  -9.30  i  0.66 

Overall  average  -9.63  i  0.40 


An  estimation  of  the  heat  of  formation  at  298® K  can  be  made  on  the  basis  of  the  assumption  that 
the  change  in  heat  content  of  the  Na-Sn  is  equal  to  the  change  in  heat  content  of  the  tin  plus  0.8  times 
the  change  in  heat  content  of  the  sodium.  Values  for  heat  capacities  are  taken  from  K.  K.  Kelley.*^ 
The  /i  H  so  calculated  is  -tO.6  kcal/mole. 
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DISCUSSION  OF  THE  WORK  ON  Na-Sn 

Although  a  conclusive  proof  of  the  method  of  analysis  cannot  be  claimed,  the  fact  that  an  average 
result  from  the  unstirred  run  is  not  far  different  from  the  average  result  from  stirred  runs,  and  the 
fact  that  the  final  average  of  the  two  sets  of  additions  includes  two  of  three  individual  results  of  each 
type  in  its  probable  error  range,  indicate  that  the  method  has  merit.  This  number  of  coincidences  is 
quite  improbable. 

A  second  indication  is  that  the  estimated  heat  of  formation  for  the  temperatures  at  which  other 
investigators  have  quoted  data  is  in  the  same  range.  Table  6  lists  the  results  of  the  other  investigators. 
The  value  of  -16  kilocalories  given  by  Biltz^  is  the  result  of  a  recalculation  of  old  data  using  corrected 
heats  for  some  of  the  intermediate  reactions.  One  reaction  in  particular, 

FeClj  +  1/2  Hg  FeClg  +  HCl  +  4H, 

was  changed  by  a  very  large  factor.  In  1924,  Biltz  and  Holverscheit^  used  the  value  10.3  kilocalories 
for  the  above  reaction  in  6.3N  HCl,  and  calculated  -10.7  for  the  heat  of  formation  of  Na-Sn.  In  1928, 
Biltz  and  Meyer  ^  recalculated  the  old  data  using  a  value  of  16.0  kilocalories  for  the  above  reaction. 

They  reported  —16.0  kilocalories  for  the  corrected  value.  No  late  data  are  available  for  this  reaction, 
so  that  a  recalculation  is  not  possible. 


Table  6.  Summary  of  heat  of  formation  data  for  Na-Sn. 


Heat 

form. 

kcal/mole 

Prob. 

error 

kcal 

Method 

Temp. 

°K 

Observer 

-11 

0.5* 

Diff. 

soFn. 

298 

Biltz  and 

Holverscheit^ 

-12.0 

0.7t 

Cal. 

298 

Kubaschewski® 
Kubaschewski  and  Seith 

CD 

1 

Diff. 

soFn. 

298 

Biltz, ^  and  Biltz  and 
Meyer® 

-0.63 

0.40 

Cal. 

873 

This  research 

-10.6t 

0.6 

Cal. 

298 

This  research 

In  the  solution  methods,  small  errors  in  either  of  the  heats  of  solution  whose  difference  is  the  de¬ 
sired  value  cause  a  relatively  large  error  in  the  result.  In  Kubaschewski’s®  experiments  a  source  of 
error  is  that  the  paraffin  coating  the  pure  sodium  would  tend  to  evaporate  at  the  high  temperature  of 
molten  Na-Sn.  This  would  leave  the  highly  reactive  surface  of  the  melt  unprotected.  Any  oxidation  of 
this  surface  would  directly  affect  the  result,  tending  to  increase  the  heat  of  formation. 

Kubaschewski®  describes  a  source  of  error  in  his  procedure,  namely,  that  the  heat  content  of  the 
crucible  above  the  temperature  of  the  calorimeter  is  only  approximately  known.  In  addition,  the  hot 
crucible  will  lose  an  almost  indeterminate  amount  of  heat  while  it  is  outside  the  calorimeter. 

The  results  of  the  present  research  would  tend  to  be  low,  as  all  of  the  approximations  made  in  the 
development  result  in  the  neglect  of  small  factors  whose  inclusion  would  tend  to  increase  the  measured 
heat  effect.  Any  increase  would  be  small,  probably  not  over  3%  of  the  calculated  value  or  about  300 
calories. 

No  definite  conclusions  regarding  the  accuracy  of  the  result  obtained  by  comparison  with  other 
values  is  possible  except  that  the  results  obtained  in  this  research  tend  to  be  lower  than  the  other  re¬ 
sults  quoted. 


♦Estimated  from  data  given  by  Biltz  and  Holverscheit. 
tCalculated  from  value  of  6%  given  by  Kubaschewski  and  Seith. 
lEstimated  value  at  this  temperature,  see  Table  5  and  following  paragraph. 
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Runs  were  first  made  with  no  stirrer  as  it  was  thought  that  the  stirring  caused  by  the  electro¬ 
magnetic  field  would  be  sufficient  to  keep  the  melt  mixed.  After  several  runs  were  completed,  it  was 
evident  from  the  reaction  rate  that  stirring  would  be  necessary.  It  was  also  hoped  that  stirring  would 
minimize  vaporization  losses  as  these  were  so  large  in  run  3  and  run  8*  that  the  runs  could  not  be 
used.  (A  contributory  fact  to  the  high  vaporization  losses  was  that  the  temperatures  in  run  3  and  run 
8  were  higher  than  in  run  1.) 

That  stirring  did  increase  the  reaction  rate  is  evident  when  results  of  run  10  and  run  11,  Table  7, 
are  compared  with  runs  1,  3,  and  8;  in  addition,  vaporization  was  reduced  markedly. 


Table  7.  Amounts  of  sodium  vaporized  in  various  runs. 


Run 

Approx, 
av.  temp.  °C 

%  Na  vaporized 

Av.  time  for 
addition,  Min. 

1 

615 

2.42 

24 

3 

660 

22.4 

22 

8 

630 

11.6 

18 

10 

660 

0.4 

11 

11 

640 

0 

13 

One  aspect  of  the  stirred  runs  that  was  not  anticipated  was  that  local  heating  was  apparently  pre¬ 
sent.  The  location  of  the  stirrer  just  below  the  thermocouple  well  is  presumed  to  have  caused  a  flow 
of  reacting  melt  to  pass  over  the  thermocouple  and  raise  its  temperature  above  the  actual  average 
temperature  of  the  melt.  No  attempts  to  avoid  this  type  of  irregularity  have  been  made,  as  the  dimen¬ 
sions  of  the  present  calorimeter  are  such  that  any  other  configuration  of  apparatus  would  necessitate 
major  changes  in  design. 

The  dispenser  was  quite  erratic  in  its  operation.  This  was  attributed  to  the  fact  that  the  sodium 
pellets  were  so  light  that  the  plates  did  not  always  operate  properly.  In  addition  the  sodium  metal  is 
quite  soft,  and  almost  any  small  force  is  sufficient  to  warp  the  pellets  out  of  shape  so  they  would  wedge 
themselves  in  the  constrictions  in  the  dispenser.  This  wedging  caused  several  runs  to  be  worthless. 

A  minor  difficulty  encountered  was  that  the  converter  was,  at  times,  difficult  to  control,  and  seem¬ 
ed  to  operate  in  a  stable  manner  at  only  a  few  power  settings.  This  made  it  quite  difficult  to  operate 
at  any  temperature  other  than  the  equilibrium  temperature  for  the  stable  settings.  The  spark-gap 
was  cleaned  twice  and  appeared  to  be  in  good  condition  each  time,  so  that  the  peculiarities  were  at¬ 
tributed  to  a  worn  control  mechanism. 

CONCLUSION  (Na-Sn) 

A  new  method  of  analysis  of  the  temperature -time  curve  has  been  devised.  Various  types  of 
curves  have  been  obtained  for  the  same  reaction,  and  the  application  of  the  method  has  given  results 
which  are  consistent  within  the  limits  of  probable  error. 

The  heat  of  formation  of  Na-Sn  was  measured  and  found  to  be  -9.63  t  0.40  kilocalories  per  mole 
for  the  reaction  at  873°  K.  On  the  basis  of  an  assumption,  the  heat  of  reaction  was  extrapolated  to  298°  K 
at  which  temperature  it  is  estimated  to  be  -10.6  t  0.6  kilocalories  per  mole. 


’•’Comparison  of  runs  1,  3,  and  8,  shows  the  following: 


Run 

Av.  T.  “C 

%  Na  vaporized 

1 

615 

2.42 

3 

650 

22.4 

8 

625 

11.6 
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PART  II 

INTRODUCTION  TO  THE  WORK  ON  Li-Sn 

As  brought  out  in  the  introduction  to  the  work  on  the  determination  of  the  heat  of  formation  of  Na- 
Sn  alloys,  the  high-frequency  induction  coil  was  first  used  as  a  calorimeter  by  Chip  man  and  GranU  in 
about  1943.  LeRoy  A.  Bromley*  patterned  such  a  calorimeter  along  the  same  lines  in  1946,  including 
suitable  control  features  which  are  described  in  the  section  on  “Description  of  Apparatus  and  Oper¬ 
ating  Procedure, Part  I. 

Where  similar  experimental  methods  have  been  employed,  the  reference  data  now  available  largely 
fail  to  present  the  equations  and  theoretical  bases  whereby  the  final  results  tabulated  were  obtained. 

In  the  work  on  Na-Sn,  a  mathematical  method  of  analysis  is  derived  enabling  one  to  interpret  the  data 
collected  using  the  high-frequency  calorimeter.  The  analytical  procedure  evolved  and  the  basic  as¬ 
sumptions  stated  will  be  found  in  the  “Derivation  of  New  Method  of  Analysis  of  Temperature -Time 
Curves,^’  Part  I,  and  other  sections  of  this  report.  This  derivation  was  used  in  calculating  the  final 
results  tabulated  in  the  summary. 

The  method  does  possess  definite  shortcomings.  This  can  be  seen  from  the  sample  calculations 
included  for  the  Na-Sn  determinations.  The  calculations  themselves  are  long  and  laborious,  especially 
so  until  one  becomes  familiar  with  the  method  of  attack.  On  the  other  hand,  shortening  the  method 
would,  of  necessity,  cut  down  one's  accuracy  of  analysis. 

Kubaschewski  and  SeithP  list  a  value  for  the  heat  of  formation  of  Li-Sn,  and  in  the  absence  of  any 
further  data,  values  determined  using  the  high-frequency  calorimeter  will  be  compared  with  theirs. 


The  real  test  of  the  method  lies  in  the  consistency  of  results  obtainable.  Barring  experimental 
error,  if  consistent  values  can  be  obtained  from  different  runs  in  which  the  reaction  exhibited  varying 
rates  of  reaction,  then  a  real  contribution  has  been  made  in  providing  a  means  of  analysis  for  such 
measurements.  This  research  endeavors  to  put  the  equations  to  a  further  test. 

Table  9.  Results  of  determinations,  Li-Sn. 

Run  Addition 

Kcal/mole  at  850“  K 

6  1 

-17.83  i  2.03 

2 

-17.36  i  2.03 

3 

-18.21  i  2.03 

-17.80  t  1.17  Average 

9  1 

-17.09  t  2.03 

2 

-17.14  +  2.03 

3 

-17.75  i  2.03 

-17.33  t  1.17  Average 

-17.57  t  0.83  Overall  average 

As  shown  in  Table  9,  the  three  additions  for  run  6  and  the  three  for  rim  9  were  averaged.  The 
final  result  tabulated  is  the  average  of  these  two  averages. 


♦Unpublished  work  of  LeRoy  A  Bromley,  Ph.D.,  instructor  in  chemical  engineering,  University 
of  California. 
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TABULATION  OF  RESULTS  ON  Li-Sn 


Table  8.  S\mimary  of  runs  and  results. 


Run 

No.  of  Li 
additions 

Heat  of  reaction 
kcal/mole* 

Comments 

1 

2 

Initial  height  of  melt  in  cru¬ 
cible  too  low.  Temperatures 
obtained  were  not  dependable. 
Stirrer  was  used. 

2 

2 

Frothing  occurred  when  Li  was 
added.  Dispenser  failed. 

Stirrer  was  used. 

3 

2 

Stirrer  stopped  when  Mo  was 
added,  changing  the  steady - 
state  temperature  of  the  system. 

4 

2 

Lithium  pellet  fell  on  conical 
shield  and  failed  to  get  into 
melt.  Stirrer  was  used. 

5 

3 

A  H  values  calculated  were  not 
consistent.!  Stirrer  was  not 
used. 

6 

3 

1st  -17.83 

Instantaneous  reaction. 

2nd  -  17.36 

3rd -18.21 

Noninstantanebus  reaction  (The 
stirrer  was  not  used  in  rim  6). 

7 

3 

Power  could  not  be  maintained 
steady.  Oxidation  of  reactants 
occurred.  Stirrer  was  not  used. 

8 

3 

- 

Oxidation  of  reactants  occurred. 
Stirrer  was  not  used. 

9 

3 

1st  -17.09 

Instantaneous  reaction. 

2nd  -  17.14 

Instantaneous  reaction. 

3rd -17.75 

Noninstantaneous  reaction  (The 
stirrer  was  not  used  in  rim  9) . 

♦Valties  corrected  to  850®  K. 
tSee  ‘  'Discussion,  ’  ’  Part  IL 


"One  probable  error  calculation  was  made.  This  calculation  was  based  on  addition  of  run  6,  which 
was  a  noiiinstantaneous  addition.  This  probable  error  was  assigned  to  the  calculated  heats  of  reaction 
for  all  the  additions  made.  A  summary  of  the  factors  considered  in  such  a  calculation  and  their  rel¬ 
ative  magnitudes  can  be  found  in  the  work  bn  Na-Sn. 
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It  is  of  interest  to  speculate  as  to  the  heat  of  reaction  at  298°  K.  Assuming  that  the  change  in  heat 
content  of  Li-Sn  is  equal  to  the  change  in  heat  content  of  tin  plus  eight -tenths  of  the  change  in  heat 
content  of  the  lithium,  one  calculates: 

-4H298  =  17.57  +  0.2  (5.23)^  =  18.62  kcal/mole 

The  total  probable  error  of  this  value  is  of  the  order  of  magnitude  of  1.5  kcal. 


DISCUSSION  OF  THE  WORK  ON  Li-Sn 

It  was  stated  in  the  introduction  to  this  work  on  Li-Sn  that  the  purpose  of  the  research  was  to  put 
the  mathematical  means  of  analysis  previously  derived  to  a  real  test.  The  fact  that  nine  runs  were 
made  before  two  good  ones  were  obtained  is  sufficient  evidence  that  there  are  many  obstacles  to  be 
surmounted  before  one  can  obtain  good  data  worthy  of  mathematical  analysis.  Extensive  preliminary 
precautions  must  be  taken  to  insure  good  runs,  and,  even  so,  experimental  errors  oftentimes  occur. 

Heats  of  oxidation  of  the  reactants  are  very  large  (e.g.,  150  kcal/mole)  compared  to  the  heats  one 
expects  in  these  runs,  and  hence  one  must  strive  to  remove  all  traces  of  oxygen  during  an  experiment. 
Runs  7  and  8  were  accompanied  by  oxidation,  as  evidenced  by  the  frothing  of  the  melt  following  lithium 
additions.  Deterioration  of  the  upper  portion  of  the  reaction  crucible  and  its  surrounding  shield  pro¬ 
vided  confirmation.  Such  oxidations  can  be  partially  avoided  by  cleaning  all  metal  and  other  surfaces 
to  remove  oxide  coatings  prior  to  the  time  the  apparatus  is  assembled.  Following  the  assembly,  the 
reaction  crucible  can  be  brought  to  red  heat  while  the  calorimeter  is  under  vacuum.  Adsorbed  oxygen 
and  moisture  are  thus  removed,  the  little  remaining  being  in  an  unreactive  state.  Argon  can  then  be 
bled  in,  the  vacuum  removed,  and  the  experiment  continued  as  described  under  the  operating  procedure. 

In  the  final  runs,  the  stirrer  was  not  employed.  Although  a  more  uniform  temperature  throughout 
the  melt  might  be  expected  with  the  use  of  a  stirrer  (e.g.,  localized  heating  partially  eliminated),  the 
geometry  of  the  apparatus  was  such  that  best  results  were  obtained  when  the  stirrer  was  not  used. 

The  comments  listed  in  Table  8  will  bear  this  out.  To  obtain  reliable  temperatures  in  the  reaction 
crucible,  a  sufficient  height  of  melt  must  cover  the  thermocouple.  Oq  the  other  hand,  the  quantities 
of  reactants  must  be  such  that  the  A  T  following  a  lithium  addition  must  not  exceed  the  limits  set  in 
the  derivation.  Using  the  calorimeter  and  dispenser  employed  previously,  certain  restrictions  were 
placed  on  the  size  of  reaction  crucible,  as  can  be  seen  from  the  above  considerations.  The  resulting 
size  actually  required  was  such  that  the  stirrer  was  stopped  following  a  Mo  addition  because  of  insuf¬ 
ficient  clearance.  This  produced  a  new  steady -state  temperature  and  made  the  entire  run  incapable 
of  calculation. 

In  run  5,  three  additions  of  lithium  were  made  and  calculated.  The  first  two  additions  gave  quite 
low  results,  and  the  last  addition  gave  a  very  high  result  compared  to  the  final  average  value  tabulated. 
Such  calculated  values  lead  one  to  surmise  that  a  fraction  of  the  lithium  pellets  were  held  up  in  the 
first  two  additions,  and  that  the  ''hold-ups''  all  went  in  with  the  last  addition.  Such  an  occurrence  is 
quite  possible. 

In  general,  the  reaction  rates  in  this  research  were  more  rapid  than  in  the  case  of  the  Na-Sn  prep¬ 
arations.  Also  the  dispenser  operation  was  much  better.  One  would  expect  this  to  be  the  case  since 
lithium,  though  less  dense,  has  a  higher  melting  point  and  is  less  subject  to  deformation  at  low  tem¬ 
peratures  than  is  sodium.  Consequently,  pellets  are  not  as  prone  to  get  "hung-up"  in  the  compart¬ 
ments  of  the  dispenser. 

Since  some  of  the  calculations  for  a  given  addition  are  subject  to  varying  interpretation,  pains 
were  taken  to  put  the  probable  error  on  the  "safe"  side.  From  the  tabulated  final  results  in  Table 
9  it  is  evident  that  such  is  the  case  since  the  actual  results  calculated  are  better  than  the  listed 
probable  error  justifies. 
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Of  the  six  additions  listed  in  the  final  tabulation  of  results,  there  is  no  discernible  trend  of  values 
for  either  instantaneous  or  noninstantaneous  type  reactions.  This  in  itself  lends  strong  support  to  the 
mathematical  procedure  derived,  which  appears  to  work  equally  well  for  either  reaction  type. 


CONCLUSION  (Li-^n) 

The  new  method  of  analysis  of  temperature -time  curves  has  been  given  another  test.  The  results 
for  both  instantaneous  and  noninstantaneous  type  reactions  are  well  within  the  probable  error  assigned. 

The  heat  of  formation  of  Li-Sn  was  measured  and  found  to  be  17.57  1  0.83  kilocalories  per  mole 
for  the  reaction  at  850®  K.  The  heat  of  formation  was  estimated  as  18.62  kilocalories  per  mole  at 
298° K.  The  probable  error  is  of  the  order  of  1.5  kilocalories  at  this  temperature.  Kubaschewski  and 
Seith®  give  a  value  of  16.8  1 1.0  kilocalories  at  this  same  temperature.  No  other  reference  values  are 
listed. 

The  probable  error  of  this  research  is  sufficiently  large  to  make  the  reference  value  and  the  val¬ 
ue  obtained  from  this  work  consistent  if  the  former  value  is  low,  and  the  latter  is  high. 


Symbol 

A 

a 

B 


B' 

C 


D 

d 

F 

f 

f'  plot 

g 

h 

*^c 

AH 

k 

L 


DEFINITION  OF  SYMBOLS 

Definition 

Area,  cm^ 

A  point  on  the  curve  of  Figure  6. 

Defined  as  radiation  heat  transfer  coefficient.  Aj  is  used  as  the  first  approximation 
to  that  quantity. 

Corrected  value  of  B. 

Integration  constant;  point  on  curve  of  Figure  6. 

Heat  capacity  of  system,  cal/^’K. 

Final  heat  capacity  of  system,  cal/  ®K. 

General  heat  capacity,  cal/®K  mole. 

Convection  heat  transfer  coefficient,  hA,(kA/L),  cal/min^K. 

Prefix  denoting  a  differential. 

Factor  defined  by:  Cp^j^y  =  Cpg^^  +  FCpj^^. 

Net  heat  transfer  rate,  -(qu  ~  qjt),  cal/min. 

The  plot  on  which  the  heat  Qjj  is  found. 

Acceleration  of  gravity. 

Convection  heat  transfer  coefficient,  cal/^C  min  cm2. 

Natural  convection  heat  transfer  coefficient;  in  McAdams  formulae,  Btu/hr  ft^  °F . 
Heat  of  Reaction;  heat  of  formation,  cal/ mole  or  kcal/mole. 

Conduction  heat  transfer  coefficient,  cai/°C  cm^  min/cm. 

Length,  feet  in  McAdams  formulae. 
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Symbol 

In 

P 

P' 

Q 

Ql 

Qn 

q 

qit 

qiT 

R 

R' 

T 

Tr 

Tl 

T2 

t 

AT 

AT' 

at" 

(AT')o 

Wt) 

w 

/3 

€ 

M 

P  ' 


Definition 

Natural  logarithm. 

Defined  as  (4B(t'  +  D);  used  as  first  approximation  to  this  quantity. 

Corrected  value  of  P. 

Heat;  or  heat  required  to  raise  alkali  from  the  dispenser  temperature  to  the 
temperature  t.,  cal. 

Heat  equivalent  to  the  lithium  vaporized,  cal. 

Heat  liberated  by  part  of  the  noninstantaneous  reaction  that  can  be  transformed  to 
ideal  instantaneous  reaction  and  is  equal  to  Cp(AT),  cal. 

Heat  liberated  by  reaction  after  ‘instantaneous”  portion.  Found  from  area  on  f  ^ 
plot  out  to  practical  time  limit,  cal. 

Heat  liberated  after  practical  time  limit  of  Qjj.  This  quantity  equals  Cp(AT). 

Rate  of  heat  transfer,  cal/min. 

Rate  of  heat  loss  at  temperature  t,  cal/min. 

Rate  of  heat  loss  at  temperature  T,  cal/min. 

Defined  as  {6B(t^)^);  used  as  first  approximation  to  this  quantity. 

Corrected  value  of  R. 

Any  temperature,  ®C,  ®K. 

Room  temperature,  ®K. 

Temperature  of  inside  turn  on  molybdenum  shield,  °K. 

Temperature  of  outside  turn  on  molybdenum  shield,  °K. 

Steady-state  temperature  before  addition,  °C,  °K. 

Equilibrium  temperature  after  addition,  °C,  ®K. 

Temperature  difference  which  would  have  been  produced  if  the  reaction  had  been 
instantaneous;  or,  in  the  Appendices,  any  temperature  difference,  ®C,  ®K. 

T  “t',  °C,  “K. 

AT ' 

Smoothed  temperature  differences  obtained  from  the  In - - - vs.  6  plots,  *C,  ‘K. 

Temperature  difference,  T  -t'  at  the  time  corresponding  to  0q,  °K, 

Temperature  difference,  T  -t'  at  the  time  corresponding  to  0,  °K. 

'  ^0  actual  curve  ^  ^0  theoretical  curve'  * 

Weight  in  moles. 

Coefficient  of  volumetric  expansion, 

Emissivity  of  a  given  surface,  no  dimensions. 

Viscosity,  in  McAdams  general  equation,  Ib/hr/ft. 

o 

Density;  in  McAdams  general  equation,  Ib/ft. 
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Symbol  Definition 

a  Stefan -Boltzmann  radiation  constant. 

e  Time,  minutes. 

0.  Time  for  the  instantaneous  reaction,  minutes. 

0Q  Reference  time  in  the  calculation  of  Qjjj. 
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APPENDIX  I 


EVALUATION  OF  D 

The  for  the  system  as  used  in  equation  (7)  can  be  represented  by  an  equation  such  as  this: 

D  =  XhA  +  XkA/L. 

Thus  evaluation  of  all  the  contributing  factors  allows  estimation  of  D. 

Assume  that : 

1.  Except  for  radiation,  natural  convection  is  the  principal  factor  for  heat  loss  from  the  sides  and 
top  of  the  crucible. 

2.  An  argon  atmosphere  is  sufficiently  like  one  of  air  to  justify  the  use  of  empirical  equations  for 
heat  transfer  in  air  to  calculate  natural  convection  losses. 

3.  Heat  losses  through  the  thermocouple  well  and  the  stirrer  are  by  conduction  only. 

4.  Heat  losses  through  the  zircon  stand  are  by  conduction  only. 

5.  Natural  convection  heat  losses  are  small  compared  to  radiation  losses,  so  that  temperatures  of 
the  outside  shield  and  base  of  the  zircon  stand  can  be  estimated  using  radiation  alone. 

The  calculation  of  D,  using  the  above  assumptions,  follows.  Calculation  of  temperature  of  the  out¬ 
side  shield: 

The  emissivity,  e',  for  an  oxidized  molybdenum  surface  is  0.7. 

aeh{t^  -Tj)  =a€A(T^  -  t|) 

<7£A(t|-t|)  =  a€'A{T4-T4) 

Calculating  e, 

€  _  _ 1 _ 

(l/£')  +(Ai/A2)  [(!/€')-!] 

°  (1/0.7)  +  [(1/0.7)  -  1]  " 

Solving,  using  T^,  as  300° K,  t  as  875° K,  and  cancelling  areas,  780° K,  and  T2  =  630° K.  Using  the 

latter  temperature,  and  an  empirical  equation  of  McAdams,® 

h^.=  0.28(4T/L)0-25, 

c 

we  find:  AT  =  (630  -  300)  =  330‘’C.  =  (9/5)(330)  =  594°F;  L  =  8.8/(12)(2.54)  =  0.289  ft;  and  =  0.28 
(594/0.289)0-25  or  h^  =  1.89  Btu/hr  ft^  °F.  This  value  is  for  natural  convection  from  the  sides  of  the 
outside  shield,  and,  expressed  in  our  units,  is  0.0154  g  cal/min  cm^  °C. 

Calculation  of  Natural  Convection  from  the  Top  of  the  Crucible 

Assume  this  to  be  equivalent  to  a  horizontal  plate  at  the  temperature  of  the  melt.  McAdams^®  gives 
an  approximate  equation  for  air  for  this  type  of  system, 

=  0.38(AT)0-25. 

AT  =  (875  -  300)  =  575‘’C;  AT  =  (9/5)(575)  =  1035°F.  =  0.38(1035)0-25  =  2.17  Btu/hr  ft2  °F  or  0.0176 

g  cal/min  cm2 
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Calculation  of  T^,  the  Temperature  of  the  Inner  Horizontal  Face  of  the  Zircon  Stand 

The  kfor  zircon  is  0.005  g-cal/sec  cm^  °C/cm,  or  0.30  g-cal/min  cm2  °c/cm. 

Assume  conduction  through  the  zircon  wafer  (0.3  cm  thick)  to  inner  horizontal  face  and  then  heat 
is  transferred  by  radiation  from  this  face  to  the  copper  stand.  Take  the  temperature  of  the  top  of  wafer 
As  875®  K,  then, 

q  =  (kA/L)(875  -T,;)  =<7€'A(T^  -  300'^); 

L  =  0.3  cm.;  e-  0.7;  substituting  and  solving,  Tj  =  840® K. 

Calculation  of  Conduction  Across  Gap  inside  Stand 

k  =  0.035  Btu/hr  ft^  °F/ft  =  0.035(60){0.00413)  =  0.0087  g-cal/min  cm2  °c/cm.  Using  the  value 
for  kj  kA/L  =  (0.3)(0.3x7rx5)/2i8,  or  kA/L  5=  0.504  g-cai/min  °C.  Since  there  is  a  very  poor  thermal 
contact  where  the  stand  meets  the  copper  block,  it  is  estimated  that  half  of  the  above  value  for  kA/L  is 
effective,  or  kA/L  =  0.252  g  cal/min  ®C. 

Calculation  of  Loss  Out  Thermocouple  Well 

The  cross-sectional  area  of  the  molybdenum  well  is  0.95  cm2,  ^he  length  to  the  junction  with  the 
copper  tube  is  15  cm.  The  area  of  the  copper  tube  is  0.30  cm^,  and  its  length  is  25  cm.  The  tempera¬ 
ture  of  the  end  Of  the  thermocouple  in  the  melt  is  875®  K,  the  temperature  of  the  end  of  the  copper  tube 
is  300® K.  The  k  for  molybdenum  is  19  g-cal/min  cm2  «c/cm,  and  the  kfor  copper  is  50  g-cal/min  cm^ 
®C/cm.  Now  the  ratio  of  the  AT's  of  molybdenum  to  copper  is: 

^'^Mo/ ^^Cu  =  ^u  ^Cu  ^Mo/^Mo  ^Mo  ^Cu  ""  50(0.30)15/19(0.95)25  =  0.50 

The  overall  AT  is  (875  -  300)  =  575°C.  We  want  the  equivalent  length  of  molybdenum  which  would  have 
a  4T  of  575®C  under  the  same  conditions.  This  equivalent  length  is  3(15.0)  or  45  cm.  Using  this  equiv¬ 
alent  length,  kA/L  =  19(0,95) /45  =  0.401  g-cal/min  °C. 

Calculation  of  the  Loss  Out  the  Stirrer 

As  an  approximate  value,  use  an  equivalent  length  of  molybdenum  shaft  of  45  cm;  k  is  19  g-cal/min 
cm^  °C/cm;  and  A  is  0.495  cm^.  Then  kA/L  is  19(0.495)/45  =  0.210  g-cal/min  ®C. 

Now  we  have  the  quantities  necessary  to  evaluate  D.  For  the  runs  without  the  stirrer,  D  =  (0.0154) 
(110.5)  +  (0.0l76)(12.6)  +  0.044  +  0.252  +  0.401  =  1.702  +  0.222  +  0.044  +  0.252  +  0.401  =  2.621  g-cal/min 
®C. 

For  rtins  with  the  stirrer,  D  =  2.62  +  0.21  =  2.83  g-cal/min  ®C. 

Check  of  assumption  (2),  McAdams  gives  an  expression  from  which  an  h^,  can  be  found  for  any  gas¬ 
eous  medium, 

For  air,  Cpfi/k  is  0.74;  for  argon,  0.69.  For  air,  p  is  ce  29;  for  argon,  a  40.  For  air,  0  is  1/T  ;  for 
argon  it  is  1/^  also.  For  air,  m  is  0.034;  for  argon,  0.042.  Then  the  ratio  of  the  above  equation  for  air 
to  argon  is,  assuming  k’s  are  equal, 

.air_  ^  (0.74)(29)^(0.042)^  ^  (1.07)(0.725)2(1.233)2  =  (0.86)0'25  =  0.963 
argon  (o.69)(40)2(0.034)2 

The  value  of  0.963  for  this  ratio  justifies  the  assumption  made,  as  the  approximate  equations  are 
probably  no  more  exact  than  0.96. 
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EVALUATION  OF  F 

Comparison  of  the  heat  capacities  of  binary  alloys  with  their  elements  indicates  that  the  molal  heat 
capacity  of  the  alloy  is  usually  less  than  the  total  molal  heat  capacities  of  its  constituents.  Using  an 
empirical  expression  for  the  heat  capacity  of  an  alloy, 

Cr,  =  Cr,  +  FC„  , 

^alloy  ^noble  metal  Phase  metal 

and  solving  for  the  factor  F,  one  ordinarily  finds  values  ranging  from  0.75  to  0.9.  Since,  in  this  research, 
a  value  for  F  is  necessary  to  enable  one  to  estimate  the  Cp  of  the  system  at  various  points  during  a 
run  a  value  of  0.80  was  chosen.  Although  this  value  is  below  the  average  of  the  calculated  values>  one 
feels  justified  in  its  use,  at  least  qualitatively,  as  sodium  is  more  electropositive  than  any  of  the  other 
base  metals  investigated.  Furthermore,  since  the  proper  value  for  F  in  this  case  is  doubtful,  a  large 
probable  error  is  assigned.*  Thus,  a  compensation  for  this  value  of  F  is  adequate. 


*The  value  of  F  used  in  computing  probable  error  is  0.8  ±0.1. 
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APPENDIX  III 

SUMMARY  OF  CALCULATION  OF  PROBABLE  ERRORS  (Na-Sn) 


Value  oiZ  Cp* 

Amount  of  sodium  vaporized  t 
Value  of  F  + 

Value  of  Cp  at  beginning  of  addition  2 

§ 

Value  of  aT  for  addition  2 
Value  of  Qi  (Cp(AT)) 

Vajue  of  Qq** 

Value  of 
Value  of 

(Value  of  Qq  +  Ql  ^11  ^IIl) 

Value  of 
Value  of  AH*** 

Value  of 


Run  1,  Addition  2 

24.14  ±  3.18  cal/"C. 

0.46  i  0.10  g. 

0.80  t  0.10 

23.01  i  3.20  cal  /  ‘’C. 

19.8  i2.85°C. 

455  i  128  cal 
192  ±  55  cal 
624  t  61  cal 
30  1  0  cal 

5.93  1  0.74  kcal/mole 
4.81  t  0.19  kcal/mole 
“10.74  "L  0.77  kcal/mole 
-10.76  +  0.77  kcal/mole 


*This  value  was  obtained  by  taking  maximum  and  minimum  reasonable  values  for  the  slope  on  the 
In  AT'  vs.  9  plot,  and  calculating  the  possible  variation  in  ^Cp.  This  variation  was  averaged  to  give 
the  probable  error. 

i  This  value  was  obtained  from  a  consideration  of  the  errors  in  weighings  and  in  the  analysis. 

^For  a  discussion  of  evaluation  of  F,  see  Appendix  II . 

^This  value  was  calculated  by  taking  the  maximum  possible  deviations  of  the  slope  on  the  in  AT' 
vs.  6  plot  and  following  the  procedure  outlined  under  '^Calculation  of  the  Heat  of  Reaction  from  Curve 
3,  Run  1,  in  Part  I  to  give  two  extreme  values  of  AT.  The  average  deviation  of  these  AT's  from  the 
quoted  AT  was  used  as  the  probable  error. 

**This  is  calculated  from  the  probable  error  in  vaporized  sodium  and  the  probable  error  in  heat 

content  values.  , 

tt  Using  the  slopes  mentioned  in  §  above,  the  f  plot  calculation  was  carried  out.  The  probable  error 

reported  is  the  average  of  the  deviations  from  the  quoted  value. 

^The  error  in  Qm  was  estimated  to  be  negligible,  since  Qjjj  was  calculated  from  an  exact  theoreti- 

cal  expression.  ^ 

§§This  value  was  calculated  from  heat  capacity  data  given  by  K.  K.  Kelley. 

♦**This  value  represents  the  combination  of  ^  +  Qj  +  Qjj 
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EVALUATION  OF  F 

Comparison  of  the  heat  capacities  of  binary  alloys  with  their  elements  indicates  that  the  molal  heat 
capacity  of  the  alloy  is  usually  less  than  the  total  molal  heat  capacities  of  its  constituents.  Using  an 
empirical  expression  for  the  heat  capacity  of  an  alloy, 

Cn  +FC„ 

i^alloy  ^noble  metal  *^base  metal 

and  solving  for  the  factor  F,  one  ordinarily  finds  values  ranging  from  0.75  to  0.9.  Since,  in  this  research, 
a  value  for  F  is  necessary  to  enable  one  to  estimate  the  Cp  of  the  system  at  various  points  during  a 
run  a  value  of  0.80  was  chosen.  Although  this  value  is  below  the  average  of  the  calculated  values^  one 
feels  justified  in  its  use,  at  least  qualitatively,  as  sodium  is  more  electropositive  than  any  of  the  other 
base  metals  investigated.  Furthermore,  since  the  proper  value  for  F  in  this  case  is  doubtful,  a  large 
probable  error  is  assigned.*  Thus,  a  compensation  for  this  value  of  F  is  adequate. 

*The  value  of  F  used  in  computing  probable  error  is  0.8  i  0.1. 
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SUMMARY  OF  CALCULATION  OF  PROBABLE  ERRORS  (Na-Sn) 


R\in  1,  Addition  2 


Value  of  £  Cp* 

24.14  ±  3.18  cal/^C. 

Amount  of  sodium  vaporized t 

0.46  1  0.10  g. 

Value  of  F  + 

0.80  t  0.10 

Value  of  Cp  at  beginning  of  addition  2 

23.01  t  3.20  cal  /  “C. 

s 

Value  of  aT  for  addition  2 

19.8  i2.85''C. 

Value  of  Qj  (Cp(AT)) 

455  ±  128  cal 

Value  of 

192  ±  55  cal 

Value  of  Qjjtt 

624  i  61  cal 

Value  of  Qjjj  it 

30+0  cal 

(Value  of  Qq  +  Qj  +  Qjj  +  Qjjj) 

5.93  i  0.74  kcal/mole 

Value  of 

4.81  i  0.19  kcal/mole 

Value  of  AH*** 

-10.74  +  0.77  kcal/mole 

Value  of 

—10.76  +  0.77  kcal/mole 

*This  value  was  obtained  by  taking  maximum  and  minimum  reasonable  values  for  the  slope  on  the 
In  AT'  vs.  Q  plot,  and  calculating  the  possible  variation  in  XCp.  This  variation  was  averaged  to  give 
the  probable  error. 

tThis  value  was  obtained  from  a  consideration  of  the  errors  in  weighings  and  in  the  analysis. 

^For  a  discussion  of  evaluation  of  F,  see  Appendix  n. 

^This  value  was  calculated  by  taking  the  maximum  possible  deviations  of  the  slope  on  the  In  AT' 
vs.  0  plot  and  following  the  procedure  outlined  under  ‘‘Calculation  of  the  Heat  of  Reaction  from  Curve 
3,  Run  1,  in  Part  I  to  give  two  extreme  values  of  AT.  The  average  deviation  of  these  AT's  from  the 
quoted  AT  was  used  as  the  probable  error. 

♦♦This  is  calculated  from  the  probable  error  in  vaporized  sodium  and  the  probable  error  in  heat 
content  values. 

tt  Using  the  slopes  mentioned  in  §  above,  the  f  plot  calculation  was  carried  out.  The  probable  error 
reported  is  the  average  of  the  deviations  from  the  quoted  value. 

^^The  error  in  Qm  was  estimated  to  be  negligible,  since  Qjjj  was  calculated  from  an  exact  theoreti¬ 
cal  expression. 

§§This  value  was  calculated  from  heat  capacity  data  given  by  K.  K.  Kelley.*^ 

♦♦♦This  value  represents  the  combination  of  (Qq  +  Qj  +  Qjj  +  Qhj)  and  Q. 
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